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Atomically sharp oxide heterostructures exhibit a range of novel physical phenomena 
that do not occur in the parent hulk compounds’ The most prominent example is the 
appearance of highly conducting and superconducting states"’’^ at the interface between 
the band insulators LaA 103 and SrTiOs. Here we report a new emergent phenomenon 
at the LaMnOs/SrTiOs interface in which an antiferromagnetic insulator abruptly 
transforms into a magnetic state that exhibits unexpected nanoscale superparamagnetic 
dynamics. Upon increasing the thickness of LaMnOs above five unit cells, our scanning 
nanoSQUlD-on-tip microscopy shows spontaneous formation of isolated magnetic 
islands of 10 to 50 nm diameter, which display random moment reversals by thermal 
activation or in response to an in-plane magnetic field. Our charge reconstruction 
model of the polar LaMnOs/SrTiOs heterostructure describes the sharp emergence of 
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thermodynamic phase separation leading to nucleation of metallic ferromagnetic islands 
in an insulating antiferromagnetic matrix. The model further suggests that the nearhy 
superparamagnetic-ferromagnetic transition can he gate tuned, holding potential for 
applications in magnetic storage and spintronics. 

In recent years the extensive efforts to create new types of heterostructures based on 
transition metal oxides have come to fruition . The strong interactions characterizing these 
materials are an appealing feature as they provide a nurturing platform for new phases^. 
Interestingly, in the most extensively studied LaAlOs/SrTiOs (LAO/STO) system, the non¬ 
magnetic insulator parent compounds give rise to conductivity'^’^, magnetism^^^^, and 
superconductivity^’^'^“^^ which are sharply tuned by the number of LAO layers in the 
heterostructure. Motivated by the complex physics arising at the heterointerfaces of band 
insulators, we studied LaMnOs/SrTiOs (LMO/STO) heterostructures, in which one parent 
compound (LMO) exhibits a very diverse phase diagram already in the bulk . In undoped 
LMO, orbital order due to Jahn-Teller distortions of the MnOe octahedra sets in at high 
temperatures of ~750 K. Subsequently, magnetic exchange between Mn^"^ ions leads to 
formation of ‘A-type’ antiferromagnetic (AFM) phase with Neel temperature of ~140 K in 
which ferromagnetic (FM) planes are aligned antiferromagnetically^^’^'^. At high doping, bulk 
LMO undergoes a phase transition into a FM metallic state^^. Hence the interplay between 
charge, orbital, and spin degrees of freedom combined with the fact that LMO has a polar 
structure similar to LAO, offers exciting potential for new emergent phenomena in 
superlattices and heterostructures . In particular, since the existence of magnetism and its 
possible origins in LAO/STO remain controversial ’ , the prospect of realization of tunable 
magnetic order in LMO/STO is of high fundamental and technological interest. 
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Here we use a scanning probe microscope based on a nanoscale superconducting quantum 
interference device (SQUID) that resides on the apex of a quartz tip (SQUID-on-tip 
(SOT))^^’^^ for imaging of the local magnetic structure in LMO films grown epitaxially on 
Ti02-terminated (001) STO substrates. Figures la-f show images of the out-of-plane 
magnetic field Bz(x,y) in LMO films of various thicknesses acquired at a height of 100 nm 
above the sample surface at 4.2 K after zero-field cooling (ZFC). A highly inhomogeneous 
Bz{x,y) is observed in all the films with a characteristic scale of 100 to 200 nm, comparable 
to the size of our SOTs (90 to 230 nm diameter, Table Tl). Remarkably, even though the 
Bzipc, y) structures look similar, the span of the local field varies by more than three orders of 
magnitude between the samples (see color bars). Figure Ig shows the quantitative analysis of 
the rms, and peak-to-peak, (inset), values of B^(x,y) vs. the thickness N (in unit 
cells (u.c.)) of the LMO film. For N < = 5, Bz and BJ^^ remain very small. However, at 

A = 6 a discontinuous change in their values by more than an order of magnitude occurs. A 
similar behavior, although on a much larger length scale, was recently reported by Wang et 
al. and interpreted as a phase transition from AFM to FM driven by the polar electronic 
reconstruction. Interestingly, the saturation value Mg obtained from global magnetization 
measurements at elevated fields (Fig. S6) shows similar N dependence (Fig. Ig blue dots) 
further suggesting a sharp onset of FM order in LMO/STO. 

Since global measurements show that LMO/STO films have an in-plane magnetization (Fig. 
S6), we have studied Bz(x,y) while progressively increasing the in-plane magnetic field //y 
(Movies Ml and M2). For each set value of Bz(x,y) was acquired either in the 

presence of //y = or after reducing By back to zero (to minimize noise and drift), with 
similar results. Figures 2a,b show an example of two consecutive Bz(x,y) images from 
movie M2 in N = 12 u.c. sample for — 160 and 161 mT, which appear to be 
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identical. To uncover the underlying magnetization process, we numerically subtract Fig. 2a 
from Fig. 2b, revealing small isolated dipole-shaped features in the differential image 
ABzix,y) in Fig. 2c and in the corresponding Movie M3. Figure 2d presents the numerical fit 
(Supplementary S2.7) of ABz(x,y) demonstrating that the dipole-shaped features are 
magnetization reversals of isolated nanoscale islands with in-plane magnetic moment m. As 
//|i is increased, isolated islands undergo moment reversals at random locations, as presented 
in Fig. 3e and Movie M4, indicating weak inter-island interactions. 

The observed magnetic behavior is profoundly different from the expected FM state. In a FM 
film, the magnetization reversal occurs through motion of domain walls (DW) separating 
regions of different magnetization orientation, until a homogeneous state is achieved at full 
magnetization. In LMO/STO heterostructures, in contrast, no DW motion is found, and 
moreover, the microscopic field Bz(x,y) in the fully magnetized state at much larger than 
the coercive field H,. is highly nonuniform and is similar to the ZFC Bz(x,y) as shown in 
Figs. 3a-d. This behavior is a clear hallmark of well-separated single-domain FM or 
superparamagnetic (SPM) islands. Indeed, the statistical analysis of iV = 8 u.c. sample in Fig. 
3f shows that the magnetic islands have a typical moment of m = 1.5x10"^ /Tg, comparable to 
magnetic moments of common SPM nanoparticles . Assuming magnetization of 4 per 
Mn atom and a magnetic layer thickness of 2 u.c. (discussed below), the islands have a 
characteristic diameter D 19 nm, smaller than our spatial resolution of 100 nm determined 
by the SOT size and the scanning height. 

Taking the experimental histogram of Fig. 3f, we simulate the ZFC and the fully-magnetized 
state by random spatial distribution of non-overlapping islands with in-plane moments m 
oriented either randomly (Fig. SI la) or fully aligned (Fig. Sllb-d). The resulting Bz(x,y) in 
Figs. Slla-d closely describe the amplitude and the characteristic length scale of the ZFC 
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state in Fig. 3a as well as of the nonuniform Bz(x,y) in the magnetized state in Figs. 3b-d, 
emphasizing the presence of well-separated magnetic islands. 

By summing the moments m of the flipped islands as function of magnetic field, we find that 
the total in-plane moment and its saturation value derived microscopically over a 

very small area (-^8x10"^ of the 5x5 mm^ sample area) well account for the global 
magnetization of the sample (Fig. 3g). The reversal process of the individual nanoscale 
islands thus provides a quantitative description of the macroscopic magnetization behavior of 
LMO/STO heterostructures. 

Our local measurements are carried out at 7 = 4.2 K well below the blocking temperature 
which can be estimated from the temperature dependence of the global magnetization (Fig. 
S7). Therefore, the fast temporal SPM dynamics are expected to be suppressed, leading to 
hysteresis in DC magnetization measurements^^ (Figs. 3g and S6). To investigate the 
existence of slow dynamics, H\\ was abruptly raised and a sequence of Bz(x,y) images were 
acquired at 380 sec intervals in N = 6 u.c. sample at constant = 80 mT (Movie M5). 
The corresponding ABz(x,y) images (Figs. 2e,f and Movie M6) reveal the presence of a 
random moment reversal process following the fast field increase and leading to a slow 
relaxation of the overall magnetization. The rate of the moment relaxation dM/dt decays 
monotonically with time (Fig. 2g) as anticipated for tunneling or thermal activation in a tilted 
potential. Observation of this nanoscale relaxation process is a direct manifestation of the 
slow dynamics of SPM islands at 7 < 7^. 

The finding of the inhomogeneous SPM state in all the samples with N > Nc over a wide 
range from N = 6 to 200 u.c. (Fig. S9) combined with the fact that samples are atomically 
flat (Fig. S5) strongly suggests an intrinsic rather than disorder-driven mechanism. The Mn^"^ 
ions, which are responsible for the LMO magnetism, have four 3d valence electrons in three 
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low energy t 2 g and one Bg orbital^^. The spin-aligned t 2 g electrons form a ‘core spin’ ofS = 
3/2, and the fourth electron, occupying the split Cg orbital, is slaved to the core spins by 
Hund’s coupling. The ground state is an ‘A-type’ AFM Mott insulator with alternating planes 
of opposite magnetization of 4 fj.g per Mn atom^^ (Fig. 4a). 

The charge reconstruction due to the polar structure of LMO leads to electron transfer to the 
LMO side of the interface from either the Mn orbitals at the top surface or from oxygen 
vacancies ’ . The number of transferred electrons grows with N > until it asymptotically 
reaches 0.5 excess electrons per u.c. (Fig. 4e). The Hund’s coupling between the Mn core 
spins and Cg electrons impedes hopping of any excess charges between AFM aligned sites, 
while FM alignment of the core spins allows excess charges to delocalize, thus lowering their 
kinetic energy^^. This double-exchange mechanism should therefore lead to a phase transition 
from an insulating AFM to a metallic FM state, giving rise to a conducting FM layer at the 
interface at high carrier concentrations (Fig. 4b). 

Our calculation reveals, however, that a new ground state with spontaneous phase separation 
emerges (SI). In this state nanoscale FM metallic islands of segregated excess electrons with 
p = 0.17 electrons/u.c. are formed at the interface within the undoped insulating AFM matrix 
(Fig. 4c). The competition between kinetic energy. Coulomb interactions, and the core spin 
exchange interaction determines the separation and size of the islands. Self-consistent 
calculation of the electrostatic potential shows that the metallic FM islands are confined to 
the vicinity of the interface, spreading over a finite thickness of a few u.c. (Fig. 4e). For 

N = 12, the islands reach D = 20 nm and m = l.bxlO"^ (Fig. 4f), in good general 

agreement with the typical size and magnetic moment found in the experiment (Fig. 3f). 
Importantly, our model predicts that the phase separation occurs for all thicknesses N > 
of LMO and the areal fraction of the magnetic islands remains below 1 (Fig. 4g). Since the 


6 



islands are of nanometer size and well separated from each other, the system always remains 
in an insulating SPM state. Intriguingly, if the excess electrons at the LMO/STO interface 
originate from Mn orbitals at the LMO surface rather than from oxygen vacancies, a double¬ 
layer of hole-doped and electron-doped islands will form as depicted in Fig. 4d. The 
saturation magnetization in Fig. Ig was calculated assuming this double-layer scenario 
combined with bulk LMO contribution of 0.2 /i^/u.c. due to canting in the AFM phase^^’^^. 
Consistently, for N < the residual in Fig. Ig and the weak Bz(x,y) in Figs. la,b 
apparently originate from AFM domains with different canting orientations^^. 

The intriguing interplay of AFM, FM, and SPM and the controllable nanoscale phase- 
separation provide the basis for novel physics and new functionalities in LMO/STO 
heterostructures. Our calculations show that a relatively small increase in carrier 
concentration can transform the insulating SPM into a conducting FM state opening the route 
to electric gate control of magnetism and magnetoresistance, with potential applications in 
magnetic recording and spintronics. LMO is thus an important addition to the oxide family of 
heterostructure materials offering novel possibilities for engineering of multifunctional multi¬ 
layers. 
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Figure 1. ZFC magnetic structure of LMO/STO of different thickness at 4.2 K. (a-f) 1.5 
X 1.5 |Lim^ scanning SOT images of the magnetie field B^ipc, y) of six samples of different 
thicknesses N = 4to 200 u.c. at -100 nm above the surface. Note more than three orders of 
magnitude change in color scale between (a) and (f). (g) Thickness dependence of the rms 
value Bz in scanning SOT images (red), global magnetization saturation value (blue), and 
theoretically calculated (blaek) which takes into aeeount the SPM and the AFM eanting 
contributions. The inset shows the peak-to-peak field variation B^^^ in SOT images vs. N 
near = 5. (h) Sehematic seanning SOT mieroscopy setup. 
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Figure 2. Field-driven and temporal magnetization reversal dynamics in LMO/STO. 

(a,b) Two consecutive images of Bz(x,y) (see Movie M2) in an Ai = 12 u.c. sample after a 
field excursion to = 160 (a) and 161 mT (b). (c) Differential image 

obtained by direct subtraction of (a) from (b) revealing magnetization reversal events of 
isolated SPM nanoscale islands (see Movie M3). Note an order of magnitude enhanced color 
scale, (d) Numerical fit to ASz(x,y) in (c) with four in-plane oriented islands using SOT 
diameter of 104 nm, scan height of 105 nm, and Ng = 2 u.c. with resulting magnetic 
moments rrii = 1.0 x 10^ /tg (0^ = 40 nm), m 2 = 2.6 x 10“’^ (O 2 = 20 nm), m 3 = 

4.7 X 10^ jtig (D 3 = 28 nm), and = 6.8 x 10^ jig (D 4 = 33 nm). (e) AB^ix.y) image in 
N = 6 u.c. sample showing the thermally activated magnetization reversals of the islands at a 
constant = 80 mT at t = 1140 s after the field ramp (see Movie M 6 ). (f) Same as (e) at 
t = 5360 s. (g) The magnetization relaxation rate dM/dt attained by vectorial summation of 
the reversal events m in each frame of Movie M 6 . Dotted line is a guide to the eye. 
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Figure 3. Saturated magnetic state and magnetization reversal process in Af = 8 u.c. 
sample, (a-d) Bz{x,y) after zero-field eooling (a), after excursion to = -250 mT (b), 

after excursion to = 250 mT (c), and in the presence of = 250 mT (enhancing 

instrumental noise) (d). Note the strong anti-correlation between (b) and (c) indicating the full 
magnetization reversal of the SPM islands. The differences in the intensity in (b-d) are due to 
a drift in the scanning height of the SOT. (e) x-y locations of the SPM reversal events with 
the color referring to the field at which the reversal occurred (see Movie M4). (f) 

Histogram of the magnetic moments m of the SPM reversals in (e) and the corresponding 
island diameters D (top axis) using = 2 u.c. (g) Cumulative magnetic moment M vs. 
applied field (red) attained by vectorial summation of all the reversal events m after sweeping 
the field from large negative value. The dashed line shows corresponding schematic 
reconstruction of M(//). Global M(H) measurement (blue, after subtraction of bare STO 
M(H)) normalized by the total sample area relative to the imaged area. 
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Figure 4. Theoretical model results and schematics, (a) Schematics of core Mn^^ spins in 
‘A-type’ AFM configuration (black arrows) in LMO at N < N^ = 5. (b) For high electron 
concentrations at the LMO/STO interface, a metallic layer of thickness N^ with uniform FM 
order (red arrows) should be formed at the interface, (c) Theoretically derived phase 
separated state for N > N^ with SPM islands embedded in AFM matrix, (d) Schematic 
representation of hole-doped and electron-doped SPM islands in the case of electron transfer 
from Mn^"^ orbitals at the top surface of LMO. (e) Calculated excess charge transfer per unit 
area q(N) = 0.5(1 — N^/N) (blue) and the thickness of electron charge layer (red) vs. 
LMO thickness N. (f) The diameter (blue) and the magnetic moment (red) of the SPM islands 
vs. N in the phase separated state, (g) The areal fraction of the SPM islands vs. N. 
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SI Theoretical model for magnetism in LMO/STO heterostruc¬ 
ture 

51.1 Charge distribution in the heterostructure 

As discussed in the main text, LMO/STO consists of an electron-doped layer within the LMO near 
the interface and a hole doped layer at the top surface. We estimate the charge density qe (per 
2D u.c.) of doped LMO layers using q{N) — 0.5(1 — Nc/N) (Fig. 4e), where we take the critical 
thickness Nc = 5 in conformity with experiment. This simple form can be obtained in the intrinsic 
polar catastrophe scenario [23]; however, here we treat it as an empirical formula. Since our model 
is electron-hole symmetric, from here on, we only refer to the electron-doped layer. 

51.1.1 Schrodinger-Poisson calculation 

The excess charges are confined close to the surface and interface due to electrostatics. However, 
they can lower their kinetic energy by delocalizing in the z-direction. We self-consistently obtain 
the spread Yg of the electrons from the interface along the z-direction by performing a Schrodinger- 
Poisson calculation, assuming a single hole-doped layer with charge +qe per 2D u.c. as a boundary 
condition at the top surface. This gives us an estimate of the layer-resolved charge distribution 
n(/), I being the layer index, and the effective single-particle potential leff(0 confines the 

electrons near the interface. The electric field (in the z-direction) between layers I and / -h 1 is 


SI 




Interface Layer 
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Figure SI: (a) Layer-resolved electron charge distribution at the LMO/STO interface. Charges 
only spread into a few layers (< 6) of LMO. (b) The effective potential that confines the 

excess charges is shown for various LMO thicknesses N. 


£{l^l + l) = £po\ + £s + £ii^ where £po\ = is the electric field due to alternating polar LaO^ 

and Mn 02 ~ sublayers, £<^ — —I'Kqejlo? the field due to the hole-doped layer at the surface and 
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( 1 ) 


is the electric field due to the Hartree potential for the charge distribution {n{l)}. Here e 18 
m is the low temperature dielectric constant of bulk undoped LMO. The potential is 

obtained by summing over the fields from the interface to the l-th layer. The kinetic energy is 
given by T-Lq = ^ where the energy dispersion e///(k) contains the ^-direction 

hopping t and the 2D dispersion in the xy-plane, eo(k) = —2t{cos kxa -\- cos kyo) ^ —At -h ta^k‘^^ 
with k = (kx^ ky). We work with spinless Fermions, as appropriate for the double exchange model 
(see below) assuming a uniform FM phase for the doped layers. The Hamiltonian 1-Lc + Hgff is 
diagonalized starting with an initial charge distribution {n{l)} and n{l) is obtained self-consistently 
via n{l) = ^F(^A(k)) where £x{^) and ipxiO^) stre the eigenvalues and eigenfunctions, 

respectively, and np is the Fermi function. The Fermi energy is determined by the charge neutrality 
constraint ^(0 = q> The results for n{l) and are shown in Fig. SL Since the charge 


density profile decays exponentially with the number of layers, to determine the number of doped 
layers we used a cut-off of n = 0.005. Below we show that the doped layers lead to a phase-separated 
(PS) state exhibiting superparamagnetism. 


SI.2 Phase separation in LMO/STO 

As discussed in the main text, the ‘A-type’ antiferromagnetic (AFM) state of undoped LMO consists 
of ferromagnetic (FM) planes that are aligned antiferromagnetically [18-20,31]. The AFM state can 
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be described by the Hamiltonian 


= -Jf Sj • + Jaf Sj 






1,11 


( 2 ) 


where i is the position of the Mn^^ ions on a simple cubic lattice with spacing a = 0.39 nm, jl 
denotes the directions in the FM planes and z> the out-of-plane AFM direction. and are the 
core spin [S — 3/2) and Cg electron spin, respectively, coupled via Hund’s coupling Jh- We work 
with Jf = Jaf = J > 0 and in the limit Jh ^ oo. The AFM in LMO is slightly canted, leading 
to a small magnetic moment ^ 0.2 /is per u.c. due to Dzyaloshinskii-Moriya exchange [19,30]. We 
incorporate this by assuming a background magnetic moment of ^ 0.2 /is per u.c. while estimating 
the saturation magnetization of the sample. 

LMO can be doped by injecting excess Cg electrons or holes, either chemically or electrostatically, 
as in the LMO/STO heterostructure. The kinetic energy of the carriers is described by the Mouble 
exchange’ model [32] 

Hkm = -t'^cos (a}aj +h.c.). (3) 

(ij) ^ ^ 


Here Qi is the polar angle of the core spin and t is the hopping amplitude of the carriers (a^). The 
above term prefers the core spins to align ferromagnetically {di — 0j), and thereby tends to induce 
metallicity. We take t — 0.3 eV and J = O.lt [32] for our calculations. 

The competition of FM double exchange with the AFM superexchange is believed to be at 
the root of the nanoscale phase separation in doped manganites [321 ESI ED. In the PS state, the 
long-range Coulomb interaction between non-uniform excess charge distributions plays a crucial 
role in determining the typical scale of the phenomenon. Indeed, using the above model and taking 
into account the Coulomb energy cost, we find that in bulk LMO the PS state is formed for doping 
below X = Xc ~ 0.1, giving rise to metallic FM islands with size ^4 — 20 nm and an excess charge 
density ^ Xc, embedded in an undoped insulating AFM matrix. Below, we estimate the energy of 
the PS state as a function of the FM area fraction pa and the radius R of the islands in the 2D 
case of the LMO/STO heterostructure. 

Kinetic energy: In the following, we estimate the kinetic energy £’kin(^ 5 Pa) of fho electrons 
within the FM island subjected to the effective confining potential V^fr(/). The kinetic energy of the 
electrons confined within an area fR? in the xy plane is obtained from 

where n = (n^, riy)] Ux^ riy being positive integers and 6///(n) contains 2 ;-direction hopping t and 2D 
particle-in-a-box energy levels eo(n) —4t-hta^7r(n^-hn^)/i?^ for a box of linear dimension ^/7TR. 
By diagonalizing JLq -\- we obtain the kinetic energy of the electrons EkmiR^Pa) as a function 
of R and the FM fraction pa- 

Magnetic energy: The formation of FM islands, while reducing the kinetic energy, leads to 
loss of magnetic exchange energy, which essentially limits the FM area fraction pa- As shown in 
Fig. S2, there are three possible A-type AFM arrangements for the LMO/STO structure. If the 
spin configurations of Fig. [S^ and Fig. [S^ are realized, then one expects to see a large magnetic 


signal from different AFM domains in the SOT scans for odd number of LMO layers for N < Nc^ in 
contrast to our observations (Fig. 1). Also, the configuration of Fig. S2 d is highly unlikely as our 
SOT measurements find that the SPM islands have in-plane magnetic moment. Therefore, for our 
calculations, we consider the spin configuration of Fig. |S2|3. In principle, the AFM configuration in 
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Figure S2: Different possible A-type AFM arrangements in LMO. The configuration in (c) is 
consistent with our measurements. 

LMO/STO heterostructure for N < Nc could be different from the A-type AFM in the bulk, e.g. 
G-type or C-type. However, the qualitative fact that we obtain an inhomogeneous SPM state for 
all N < 200 will not change if we take G-type of C-type AFM states as FM tendencies will be even 
more suppressed. 

For N > Nc^ Nq layers get doped with electrons. If these layers host FM islands in an AFM 
matrix with a FM area fraction pa^ then the magnetic energy of the layers is given by £’mag(Pa) = 
— (3A^e ~ ‘^Pa^e — 1) • As in the case of bulk LMO, the competition between kinetic double 

exchange and magnetic superexchange gives rise to a PS state with < 1- However, as the excess 
charges segregate within the FM regions, it costs a lot of Coulomb energy to form a large FM 
region. This essentially limits the size of the FM islands. 



N (u.c.) N (u.c.) 


Figure S3: (a) Comparison of energies of the FM, and phase separated, F^ps, states as a 

function of LMO thickness showing the stability of PS state for all N > Nc. (b) The charge density 
p{N) within a FM island in the PS state. 
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Coulomb energy: To obtain the Coulomb energy cost, we approximate the hole-doped layer 
at the surface as a uniformly charged 2D plane with surface charge density ao = qe/a^ and the 
electron doped layer at the interface as a square lattice of 2D disks, with radius R and surface 
charge density af = —cro/pa, having average spacing The Coulomb energy is obtained 

from £'couiomb = where k = (k||,fc^), is the area of the system, 

and p(k) is the Fourier transform of the 3D charge density. For FM area fraction pa < the 
Coulomb energy (per 2D u.c.) contribution from the non-uniform part of the charge distribution 
is obtained as 


-^Coulomb 



( 4 ) 


where g = 2y^(^ix + ^2y), 91^92 being integers, and Ji{x) the Bessel function, and V = jea 
is determined by the dielectric constant e = eps in the PS state. Since eps is not known, we take 
for our calculation eps ~ 100, the value for doped LMO m- However, our results do not change 
qualitatively over a range of eps values. 


SI. 2.1 Results 


Summing over E^agipa), E]^i^{R,pa), and EcoulombiR,Pa), we obtain the energy Eps{pa,R) of the 
PS state and minimize it to obtain the optimal diameter D and area fraction pa of the FM islands, as 
shown in figures. 4f and 4g. The magnetic moment m (Fig. 4f) of the FM islands is obtained from 
their volume ivR^Nea assuming d/ip per Mn atom. The total magnetic moment M of the sample 
(Fig. Ig) is calculated by summing the magnetic moments m of the electron- and hole-doped layers 
over the 5x5 mm^ area of the sample, as well as the background contribution of 0.2/iB per Mn 
for the {N — 2Ne) + 2(1 — Pa)R^e undoped AFM part of the LMO layers. Energies of the SPM and 
FM states are compared in Fig. [S^ . We find the SPM state to be stabilized over uniform FM, 
i.e., Pa < 1, for all thicknesses Q < N < 200, in conformity with our SOT measurements. The 
charge density inside each FM island varies weakly with for > 6 and stays around 0.17 (Fig. 
SS d). Figure 4f shows that the size of the FM islands is on the nm scale, giving rise to the SPM 
behavior. The calculated moments and diameters of the FM islands are in good agreement with 
corresponding typical values, D 19 nm and m 1.5 x 10^ /ip, found experimentally (Fig. 3f). 
However, in reality, disorder can give rise to a distribution of these quantities, as seen in figure. 
3f. The quantities D, m, and pa show non-monotonic dependence on A^, peaking at A^ 12 (Fig. 
4f,g). Around this thickness, a transition from insulating SPM to the metallic FM state could be 
induced by increasing the carrier concentration at the interface by an external gate voltage. 


S2 Experimental details 

All the measurements were performed at 4.2 K in He exchange gas at ^ 1 mbar. The pixel size of 
all the SOT images shown in the main text is 5 x 5 nm^. Acquiring each image took ^ 5 minutes. 


S2.1 SQUID-on-tip (SOT) characteristics 


The scanning SOT microscopy technique, including the Pb SOT fabrication and characterization, 
is described in Refs. 25,26,35. Figure S4 shows the measured quantum interference pattern Ic{H±) 
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Figure S4: Quantum interference pattern of a Pb SOT. Critical current Ic{H^) of the SOT used 
for measurement of the 8 u.c. sample vs. out-of-plane magnetic field at 4.2 K. 


of the Pb SOT used to investigate the 8 u.c. sample, which is typical for our devices. It had an 
effective diameter of 114 nm (204 mT modulation period), 66 /iA critical current at zero field, and 
white flux noise (at frequencies above a few hundred Hz) of 200 A different SOT of 

^ 100 nm diameter was used for each sample to study the local H^(x,y), as summarized in Table 
Since the 4 and 5 u.c. samples produced a very weak signal, a larger SOT of 229 nm was used 


T1 


for both samples. 

SOTs are sensitive only to the out-of-plane component of the magnetic field Bz and can operate 
in the presence of elevated in-plane and out-of-plane fields. The field sensitivity of a SOT arises from 
the field dependence of its Ic{H^) and is maximal around the regions of large \dIcldH\. Therefore, 
the SOTs usually have poor sensitivity at — 0, as seen from Fig. Using a vector magnet, 
we have applied a constant to bias the SOT to a sensitive region and then imaged the local 
Bz{x^ y) at various values of up to our highest field /io^|| = 250 mT. The presence of did not 
cause any observable effect on Bz{x^y) because of the in-plane magnetization of LMO with large 
anisotropy. The values of the applied H± for the various samples are listed in Table |TT] along with 
the estimated scanning height h of the SOT above the sample surface. For 6 to 24 u.c. samples, 
we have a more accurate evaluation of /i, obtained from the best fit to AH^(x, ^), as demonstrated 
in Fig. 2d and described in section [S277 


Sample (u.c.) 

4 

5 

6 

8 

12 

24 

200 

SOT diameter (nm) 

229 

229 

101 

114 

104 

90 

111 

HoH± (mT) 

10 

20 

28 

28 

65 

142 

65 

h (nm) 

~100 

~100 

80 

105 

105 

137 

~150 


Table Tl: SOT parameters for various samples. Listed are the SOT diameters, the applied out-of- 
plane field yoH± at the working point, and the estimated height of the scanning SOT above the 
sample surface. 


S2.2 Sample fabrication and characterization 

The LMO films on Ti02-terminated single crystal STO (001) substrates were deposited by pulsed 
laser deposition (PLD) from a polycrystalline LMO target in an oxygen partial pressure of 10“^ 
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Figure S5: Surface topography of the LMO/STO samples. Atomic force microscopy of the 12 u.c. 
(a) and 24 u.c. (b) samples showing single atomic step terraces. 


mbar at 750°C by using 1.8 J/cm^ pulses at 248 nm with a repetition rate of 2 Hz. The STO 
substrates (CrysTec GmbH, Berlin) of 5x5 mm^ and 0.5 mm thickness were double-side polished 
and chemically treated in buffered hydroffuoric acid and annealed at 950°C in oxygen, resulting 
in singly-terminated STO surface with atomically ffat terraces of single STO unit-cell height and 
terrace width of ^ 300 nm. The layer-by-layer growth of the films was monitored in situ using reflec¬ 
tion high-energy electron diffraction (RHEED) and the samples were cooled to room temperature 
in oxygen at the deposition pressure. Figure presents representative atomic force microscopy 
images of = 12 and 24 u.c. LMO films showing that atomic-step terraces are preserved, demon¬ 
strating the layer-by-layer growth of the LMO thin films. 

S2.2.1 Global magnetization measurements 



Figure S6: Global magnetic moment vs. field measurements at 4 K. (a,b) The in-plane magnetic 
moment M{H) of 5 x 5 mm^ LMO/STO samples of various indicated thickness vs. the applied 
in-plane field, (c) Saturation magnetic moment of the samples Ms (after subtraction of bare STO) 
(blue) and Mg per Mn atom (red) vs. thickness N. 
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Global magnetization measurements of the samples were done using a Quantum Design magnetic 
properties measurement system (MPMS) vibrating sample magnetometer. Figures S6r,b show the 
magnetic hysteresis M{H) loops for LMO samples of different thickness N. The ‘STO’ curve 
refers to a bare STO substrate that went through the same process, not including PLD. The finite 
hysteretic signal of the bare STO may either arise from an artifact such as residual magnetic field of 
the magnetometer’s superconducting magnet [36] or from silver paint contamination of the substrate 


The N — A and 5 u.c. samples show a very small change in magnetization relative to the bare 
STO, while a substantial difference is observed upon increasing the thickness by a single u.c. to 
= 6, as shown in Fig. S6:. The saturation magnetic moment Ms (as well as the coercive field) 


X 10“® 



x10“® x10“® 






Figure S7: Temperature dependent in-plane magnetic properties of N=12 u.c. sample, (a) Field 
cooled (FC) and zero-field cooled (ZFC) in-plane magnetic moment M{T) measured in different 
applied measurement fields /io7f||. (b) Tmax vs. /ioi7|| showing decrease of Tmax with field. (c,d) 
M{H) loops at various temperatures showing the decrease in the coercive field with temperature. 
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increases monotonically with N > Nc = 5. Figure [S^ also presents Mg per Mn atom, which shows 
a sharp jump at N = 6 and a non-monotonic behavior at larger thicknesses. The magnetization 
per Mn atom is always smaller than the expected tfis indicating that only a fraction of the Mn 
atoms are in the FM state. 

The temperature dependence of the in-plane magnetic properties of = 12 u.c. sample are 
shown in Figure revealing the onset of magnetism below Tc=100K. Field cooling (FC) was done 
using a cooling field = 1 T and a measurement field of /io^|| = 0.1 T was applied during 

the warm-up process. Zero field cooling (ZFC) measurements were done during warm-up in the 
presence of the indicated measurement field values. As shown in Fig. [S7^ , ZFC curves display a 
maximum at T^y^ax which decreases with H\\ as summarized in Fig. [S7b. In addition, magnetic 


hysteresis loops (Figs. [S^ and S71) acquired at different temperatures show that the coercive field 
/ioi^c decreases with increasing temperature, down to 10 mT at 55 K. The behavior of Tmax and the 
hysteresis loops point to a possible existence of a blocking temperature Tg >80 K 


S2.3 Additional Bz{x^y) and ABz{x^y) images 

We explored several different regions of the samples with no qualitative differences. Fignre [S8| shows 
a large area Bz{x^y) scan of 10 x 10 of the A = 12 u.c. sample after ZFC, demonstrating the 
relative uniformity of the magnetic features. 

Figure S9 shows examples of the differential ABz{x^y) images in various samples. All the 
samples with A > Ac = 5 show clear dipole-like features of SPM reversal events. For the N = 200 
u.c. sample, our maximal /ioA|| = 250 mT was insufficient to reach He in order to study SPM 
reversals. 



Figure S8: Large area image of A = 12 u.c. sample. Bz{x^y) image of 10 x 10 yuP after ZFC. The 
color scale spans 2.8 mT. 
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Figure S9: Representative ABz{x^y) images in various LMO/STO samples. Dipole-like SPM re¬ 
versal features are observed in all samples with N > Nc = 5 u.c. The images were attained by 
subtracting consecutive Bz{x^y) images with 1 mT applied field intervals. 


S2.4 In-plane anisotropy 


By applying at different angles, we find a significant in-plane magnetic anisotropy of the SPM 
islands. For i7|| oriented close (0 — 7°) to the [100] STO direction (x-axis), the angular distribution 
of the SPM magnetization reversals is peaked at 0 = 0, as shown in Fig. |S10^ and illustrated by 
the ABz{x^y) image in Fig. SlOl. For i7|| at 52°, most of the events are still oriented around 
0 = 0° (Figs. S103,e). However, few events appear at angles close to 0 = 90°. When H\ 


IS 


at 97° (Figs. |S10| :,f), the angular distribution shows a broad maximum around the y-axis ([010] 
STO). The in-plane magnetization thus shows fourfold anisotropy with fourfold easy axes along 
the LMO crystallographic directions that are locked to the underlying STO crystal structure. The 
observed differences in the anisotropy barrier for the two orthogonal directions is caused apparently 
by symmetry breaking at the cubic-to-tetragonal transition of STO at T < 105 K, leading to domain 
structure 


S2.5 Simulations of Bz{x,y) 

To simulate the local field distribution Bz{x^y) in the SPM state, we use a random distribution 
of non-overlapping magnetic islands using the experimentally attained distribution of sizes and 
moments m in the 8 u.c. sample shown in Fig. 3f. Figure |S11^ shows the ZFC state in which 
the in-plane moment orientation was taken to be random along x and y easy axes, as discussed 
in the previous section. The shown Bz{x^y) was calculated using the experimental values of SOT 
diameter of 114 nm and scan height of 105 nm. The result of Fig. |S11^ compares well with the 
experimental data in Fig. 3a both in the size of characteristic features and in the span of H^(x, y). 
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Figure SIO: (a-c) Histograms of the angular distribution of the SPM moment m orientations in 
= 12 u.c. sample for three orientations of the applied field 0 = 7° (a), 52° (b), and 97° (c) 
relative to the [100] STO orientation, (d-f) Examples of corresponding ABz{x,y) images showing 
various moment orientations. 


Figure [ST^ shows Bz{x^ y) of the same ‘sample’ as in Fig. Sll r but with all the moments oriented 
randomly either in -\-x or -\-y directions in order to simulate the fully magnetized case at > He 
applied at 45° with respect to the x axis. In contrast to the FM case, in which a uniform field is 
attained at full magnetization, the resulting Bz{x^y) remains highly inhomogeneous because the 
SPM islands are well separated. Figure |S11| 3 shows that a fully magnetized SPM state displays 
Bz{x^ y) that is similar to the ZFC state with a moderate reduction in the field span consistent with 
the experimental data in Figs. 3a-d of the main text. An SPM state fully magnetized in -\-x and 
+y directions (Fig. |S11| 3) results in Bz{x^y) that is identical to Fig. Sll3 but of opposite polarity, 
as observed experimentally in Figs. 3b,c. Similar result is attained by polarizing all the islands 
in the +x direction (Fig. indicating that all the islands have fully reversed their magnetic 

moments. 
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Figure Sll: Numerical simulations of y). (a) Numerical simulation of ZFC state with random 
position and magnetization orientation of SPM islands along -\-x^-\-y easy axes using m distribution 
of Fig. 3f. (b-d) Same as (a), with all the moments oriented randomly in the directions x^y (b), 
-\-x^-\-y (c), and only along -\-x (d) representing a fully magnetized state at » He when the 
field is applied at 45°, 225° and 0° respectively. 

S2.6 Movies of Bz{x^y) and ABz{x^y) 

Movie Ml shows a sequence of 1.5 x 1.5 /im^ Bz{x^y) images in the N = 8 u.c. sample acquired 
upon increasing from 0 to 150 mT in steps of 1 mT applied at 0 = 52° relative to the x 

axis (and the [100] STO direction) after full magnetization of the sample at —250 mT. The first 
frames of the movie mainly show the instrumental x^y drift arising from application of i7||. With 
increasing subtle changes in y) at random locations begin to be visible. These changes 

grow significantly on approaching He — 95 mT followed by reduction in the changes at higher fields. 
Note that the Bz{x^ y) images at the beginning and at the end of the movie are practically inverted, 
indicating that all the SPM islands have fully reversed their moments. 

Movie M2 shows a similar process in the N = 12 u.c. sample upon increasing from 0 

to 250 mT at 0 = 7°. After the initial drift, random small changes become visible. Our maximal 
in-plane field of 250 mT is, however, insufficient to reach a full inversion of all the SPM islands. 

The corresponding ABz{x^y) images in the = 12 u.c. sample are shown in Movie M3 for 
from 125 to 250 mT obtained by subtraction of consecutive images in Movie M2 (note an 
order of magnitude smaller color bar span). Randomly-appearing dipole-like features show the 
magnetization reversal process of SPM islands. Most of the dipole-like features are oriented close 
to the X direction along the easy magnetization axis of [100] STO. 

Movie M4 shows the x,y coordinates of the SPM island reversal events (after drift correction) 
in the N = 8 u.c. sample as is increased. The events show random uncorrelated behavior 
consistent with the SPM state. The compilation of all the locations is presented in Fig. 3e. 

The thermally activated process of the magnetic reversal of the SPM islands is presented in 
Movies M5 and M6 of = 6 u.c. sample. The in-plane field was rapidly ramped to 60 mT and 
a sequence of Bz{x^y) images (Movie M5) was acquired at a constant /io7f|| = 60 mT for about 
1.5/i. The acquisition time of each image was 360 sec with 20 sec interval between the images. 
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The sequence of ABz{x^y) images (Movie M6) is attained by subtraction of consecutive Bz{x^y) 
images. The dipole-like features in Movie M6 show that following the rapid increase in /io^|| the 
islands continue to reverse their moments predominantly in the direction of /io^|| for an extended 
period of time through a thermally activated process. The number of the reversal events decays 
with time as seen in Movie M6. The relaxation rate of the total in-plane magnetic moment in 
the scanned area presented in Fig. 2g is obtained by vectorial summation of the reversing moments 
m along /io^|| direction for each frame of Movie M6. 

S2.7 Data analysis of ABz{x,y) and the fitting procedure 

For each sample, a sequence of Bz{x^y) images was acquired, increasing in steps of 1 mT 

from 0 to 250 mT after sweeping the field from -250 mT. After numerically detecting and correcting 
the relative drifts in the x^y plane, each pair of images was subtracted to obtain ABz{x^y) (Fig. 
2c). Each dipole-like feature in ABz{x^ y) is assumed to originate from a SPM island that reversed 
its in-plane magnetic moment oriented at an angle 9 from —m to +m. The value of m is determined 
by the diameter D and thickness (taken from theory in Fig. 4e) of the island using magnetization 
of 4 /i^/Mn with lattice parameter of 0.39 nm. Fitting is then performed with 0, and h as free 
parameters including a convolution with the SOT size (keeping the same and SOT height h for 
all islands). Figure 2d shows an example of the fit result. Since D is smaller than h and the SOT 
diameter, the resulting ABz{x^ y) is essentially described by a point-like in-plane magnetic moment 
m almost independent of its size D. Therefore, the choice of (taken from theory) affects the 
derived D but has little effect on the derived value of m. 


References 

[31] Cohn, J. L., Peterca, M. &; Neumeier, J.J., Low-temperature permittivity of insulating per- 
ovskite manganites, Phys. Rev. B 70, 214433-(l-6) (2004). 

[32] Dagotto, E. Nanoscale Phase Separation and Colossal Magnetoresistance (Springer-Verlag, 
Berlin, 2002). 

[33] Nagaev, E. L. New type of self-localization state of carriers in an antiferromagnetic semicon¬ 
ductor. Pis’ma Zh. Eksp. Teor. Fiz. 55, 646 (1992). 

[34] Kagan, M. Y. & Kugel, K. I. Inhomogeneous charge distributions and phase separation in 
manganites. Physics-Uspekhi 44, 553 (2001). 

[35] Finkler, A., Vasyukov, D., Segev, Y., Neeman, L., Lachman, E. O., Rappaport, M. L., Mya- 
soedov, Y., Zeldov, E., Huber, M. E. Scanning superconducting quantum interference device 
on a tip for magnetic imaging of nanoscale phenomena. Review of Scientific Instruments^ 83, 
073702 (2012) 

[36] Quantum Design 2009 Application Note 1070-207: Using PPMS Superconducting Magnets at 
Low Field. 

[37] Colmar, F., Mudarra Navarro, A. M., Rodriguez Torres, C. E., Sanchez, F. H., Saccone, F. 
D., dos Santos Claro, P. C., Benitez, G. A., Schilardi, P. L. Extrinsic origin of ferromagnetism 


S13 



in single crystalline LaAlOs substrates and oxide films Applied Physics Letters^ 92, 262503 
(2008). 

[38] Wernsdorfer, W. (2001) Classical and Quantum Magnetization Reversal Studied in Nanometer- 
Sized Particles and Clusters, in Advances in Chemical Physics, Volume 118 (eds 1. Prigogine 
and S. A. Rice), John Wiley & Sons, Inc., Hoboken, NJ, USA. doi: 10.1002/9780470141786.ch3 

[39] Shinde, S. R., Ogale, S.B., Higgins, J.S. Zheng, H., Millis, A. J., Kulkarni,V. N., Ramesh, 
R.,Greene, R. L., Venkatesan, T. Co-occurrence of Superparamagnetism and Anomalous Hall 
Effect in Highly Reduced Cobalt-Doped Rutile Ti02 Films Physical Review Letters^ 92, 166601 
(2004) 

[40] Chen, Q., Rondinone, A.J., Chakoumakos, B.C., Zhang, Z.J. Synthesis of superparamagnetic 
MgFe 204 nanoparticles by coprecipitation Journal of Magnetism and Magnetic Materials^ 194, 
1-7 (1999) 

[41] Zhang,Y.D., Budnick, J.I., Hines, W. A., Chien, C.L., Xiao, J.Q. Effect of magnetic field 
on the superparamagnetic relaxation in granular Co-Ag samples. Applied Physics Letters^ 72 
2053-2055 (1998) 

[42] Bitoh, T., Ohba, K., Takamatsu, M., Shirane, T., Chikazawa, S. Field-cooled and zero-field- 
cooled magnetization of superparamagnetic fine particles in CU 97 C 03 alloy: comparison with 
spin-glass Au 96 Fe 4 alloy. Journal of the Physical Society of Japan^ 64, 1305-1310 (1995) 

[43] Honig, M., Sulpizio, J. A., Drori, J., Joshua, A., Zeldov, E., Ilani, S. Local electrostatic imaging 
of striped domain order in LaAlOs/SrPiO^. Nature Materials 12, 1112-1118 (2013). 


S14 



